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ABSTRACT 

High redshift Type la supernovae (SNe la) are likely to be gravitationally lensed 
by dark matter haloes of galaxies in the foreground. Since SNe la have very small 
dispersion after light curve shape and colour corrections, their brightness can be used 
to measure properties of the dark matter haloes via gravitational magnification. We 
use observations of galaxies and SNe la within the Great Observatories Origins Deep 
Survey (GOODS) to measure the relation between galaxy luminosity and dark matter 
halo mass. The relation we investigate is a scaling law between velocity dispersion and 
galaxy luminosity in the i?-band: a = a^,{L/L^,)'''' , where L, = 10^°/i~^Lq. The best- 
fitting values to this relation are cr* — 136 km s~^ and 77 — 0.27. We find cr* < 190 km 
at the 95 per cent confidence level. This method provides an independent cross- 
check of measurements of dark matter halo properties from galaxy-galaxy lensing 
studies. Our results agree with the galaxy-galaxy lensing results, but have much larger 
uncertainties. The GOODS sample of SNe la is relatively small (we include 24 SNe) 
and the results therefore depend on individual SNe la. We have investigated a number 
of potential systematic effects. Light curve fitting, which affects the inferred brightness 
of the SNe la, appears to be the most important one. Results obtained using different 
light curve fitting procedures differ at the 68.3 per cent confidence level. 

Key words: gravitational lensing - supernovae: general - dark matter 



1 INTRODUCTION 

Weak gravitational lensing can lead to stretched or mag- 
nified (or de-magnified) images of distant sources. The 
stretching of an image depends on the shear of the 
lens, whereas the magnification is related to the con- 
vergence of the lens. Measurements of shear through 
the ellipticity of galaxies lensed by galaxies (galaxy- 
galaxy lensing) have successfully be en used to probe 



1984 



the nature of dark matter halo es (iTyson et al.l 
Brainerd. Blandford fc Smaill 1 19961: Hudson et al.l 
Fischer et al.1 I2OO0I: iGuzik fc SeliaM 2002; H oekstra etlo] 



1998 



2003i. '2004': ' Kleinheinrich et all 1200 6!). In this paper we use 
measurements of convergence rather than shear to investi- 
gate properties of dark matter haloes. 

The convergence can be obtained from the bright- 
ness of standard candles, because magnified standard can- 
dles are brighter than de-magnified ones. As standard 



candles we use SNe la, which after light curve shape 
and colour corrections have small dispersion. Gravita- 
tional lensing adds extra scatter to th e Hubble dia- 
gram ('Kantowski, Vaug han fc Brancbl 19951: lFriemanlll997l : 
IWamb sganss ot al. 1997:: IHoIz fc Waldl " 1 19981 ) and is con- 
sequently regarded as a systematic uncertainty in the 
context of SN la cosmology. This e xtra dispe rsion. 



which increases with reds h ift (see, e.g., Bergstrom et al.l 
I2OOOI : iHolz fc Lindeil l2005l : iGunnarsson et al.' '2006|) can 
also be used as a c osmo logical probe (Mctcalf Il999l : 
iDodelson fc VaUinottol I2OO6I ). The distribution of SN la 
Hubble diagram residuals, which is assumed to trace the 
magnification probability distributi on, can be u sed as an- 
other gravitational lensing probe l|Wanel |2005| ). which is 
sensit ive to the fraction of compact objects in the Uni- 



verse (lRauchll99ll:lMetcalf fc Silk 1999!: ISeliak fc Hold 



Mortsell. Goobar fc Bergstrom 
Metcalf fc Silkll2007h . 



2001i : iMintv et al 
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Here we use yet another method which utilises indi- 
vidual SNe la rather than the distribution of residuals. We 
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cross-correlate SN la brightnesses and gravitational magni- 
fications computed usin g proper t ies of galaxies in the fore- 
ground, as suggested by iMetcalj (|l998l ). For an alternative, 
potentially very powerful method, also utilising correlations 
between foregr ound galaxies and supernova brightnesses, see 
iMetcalj (|200lh . 

For our investigation we use a sample of SNe la which 
was obtained within the Great Observatories Origins Deep 
Survey (GOODS). Owing to the HST observations, which 
were part of the GOODS, the sa mple contain some of th e 
most distant SNe la ever observed (|Riess et al.ll2004l200'if ). 
The high redshift of the supernovae means that they are 
likely to be significantly magnified by the matter in th e fore- 
ground (|Bergstrom et al.ll2000l : IHoIz fc Lindej|2005h . This 
magnification is then used to infer the properties of the lens- 
ing galaxies. 

The outline of the paper is as follows. In section [2] we 
describe the sample of SNe la and the galaxy catalogues we 
use. The method is explained in section [S] The results ob- 
tained when applying the method to the data are presented 
in section |4] Systematic effects which could affect the results 
are investigated in section (5] Section |6] is devoted to a par- 
ticularly important systematic effect, namely the differences 
in the results which arise due to different light curve fitting 
packages. The results are finally summarised and discussed 
in section [7] 

Throughout the paper, we assume a flat universe dom- 
inated by dark matter and a cosmological constant. Unless 
otherwise stated = 0.3 and Q.a = 1 — f^M- 



2.1.1 Hubble diagram residuals 

In order to compute the Hubble diagram residuals, AmsN, 
which tell us whether a SN la is brighter or dimmer than 
average, we have to predict the average SN la brightness 
(i.e. distance modulus) as a function of redshift. To predict 
the distance modulus we use the formula 



fj.p{z) ^ M + 51og^o D-l{z; Qm), 

where the //o-independent luminosity distance. 



Dl{z;Qm) = (1 + ^) 



dz' 



y^nmil + z'y^ + il-QM) 



(2) 



(3) 



depends only on Qm, since we assume the Universe to be 
flat, and the redshift. We use the best-fitting value of the 
norm alisation constant A4 to the Gold sample (|Riess et al.l 
l2007t ) obtained for a fixed value of Om . Inclu ding also the 
SNe la classified as silver in iRiess et all l|200'iD makes a neg- 
ligible difference to M of -0.008 mag. The Hubble diagram 
residual of a SN la is thus computed by 



AmsN 



^^0 — Mpi 



(4) 



where fio is the ex tinction cor r ected distance modulus of 
the supernova from lRiess et al.l l|2007l ). We assume the un- 
certainty in the Hubble diagram residual, (TAmq M , to be the 
same as the uncertainty in the distance modulus (iRiess et al.l 
|2003). The supernova data and the residuals computed for 
Qm = 0.3 are hsted in Table [T] 



2 DATA 

In this paper we make use of data observed within the 
GOODS project. Two fields were observed within this 
project: GOODS-North {12'^ 36™ 55% +62° 14' 15") and 
GOODS-South (3'' 32"" 30=, -27° 48' 20"). The observations 
include both galaxies and supernovae. 



2.1 Supernovae 

Thanks to the HST observations, which formed part of 
the GOODS project, some of the furthest supernovae 
ever disco vered were located in GOO DS-S and GOOD S-N 
ijRiess et a l. 2004; Stro lger et al..2004l:lRiess et al.l200"7h . We 
also include SN1997fr (|Riess et al.ll200ll l at 2 ~ 1.8 located 
in GOODS-N in our sample, although it was not detected 
in the actual GOODS search. Since the effects of gravita- 
tional lensing increase with redshift, the GOODS SNe la are 
well suited to our purposes. For each SN la we use redshift, 
distance mod ulus corrected for d ust extinction, and its un- 
certai nty from lRiess et al.l (|2007l ). The SNe la in lRiess et al.l 
(|2007l ) are classified as either gold (high confidence) or sil- 
ver (likely but uncertain). We work primarily with the gold 
SNe la. There are 7 and 19 gold SNe la located in the 
GOODS-S and GOODS-N, respectively. Our Gold sample 
hence consists of 26 SNe la. In order to compute the magni- 
fication, we need the redshift, zsn, and the position on the 
sky, dsN, for all SNe la. We will refer to these variables as 



AsN = {zsn, Ssn}, 
in the following. 



(1) 



2.1.2 Field boundaries 

Since unobserved galaxies outside of the fields can contribute 
to the magnification, only SNe la located further than a dis- 
tance Oc on the sky from the boundaries of t he fie lds are 
included in the analysis. In lGunnarsson et al.1 (|2006l ) it was 
shown that it is sufficient to include galaxies within a circular 
radius of 60" when estimating the magnification of SNe la 
due to foreground galaxies. We therefore use Oc = 60" as our 
cut-off radius. HST04Sas is located too close to the bound- 
ary of GOODS-N and is therefore excluded from the Gold 
sample. HST05Koe is located in a masked region, which lack 
observations of foreground galaxies, and is for that reason 
also excluded from the Gold sample. We consequently use 
only 24 of the SNe la in the Gold sample. 



2.2 Galaxies 

To model the dark matter haloes in the foreground of 
the SNe la we use galaxy observations obtained within 
the GOODS (Giavalisco et al. 2004; Capak et al. 2004]). 
For GOODS-S we use a VLT/ISAAC selected catalog 
with a lOg point sou rce sensitivity of ttiab = 25.1 mag 
(|Giavalisco et al.ll2004l ). while we for GOODS-N use a SUB- 
ARU/Supreme cam _R-band selection with a 5cr point source 
sensitivity of rriAB = 26.6 mag Capak et al.l (^004). The 
GOODS-N and the GOODS-S have been imaged in many 
wavelength bands, ranging from the ?7-band to the iC'-band, 
which make it possible to calculate accurate photometric 
redshifts for the galaxies in these fields. 
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2.2.1 Photometric redshifts 

Since our galaxy catalogues contain several thousands of 
objects, most of them have only photometric redshifts. To 
compute the photometri c redshifts we use a version of the 
template-fitting method (|Dahlen et al.||200 ?j). This method 
gives us for each galaxy a photometric redshift distribution, 
an absolute B-band magnitude, and the best-fitting spectral 
type. The probability distribution has sometimes multiple 
peaks, reflecting different possible solutions. Depending on 
the redshift of the SN la, different peaks could place the 
galaxy in the foreground or background of the supernova. 
The photometric redshift distributions, p{z), are in most 
cases asymmetric and the mode, i.e. the highest peak, is 
usually different from the average redshift. As our photo- 
metric redshift estimate for the galaxies in the line of sight 
we generate (for each one) a r andom number f rom the de- 
rived probability distribution. IWittmanI (|2009l ) has shown 
that this is a more robust approach than using the mean or 
median of the distribution. 

In order to compute the magniflcation of a SN la we 
need properties of the A'^gai galaxies along the line of sight. 
These properties will in the following be referred to as 



■^los — {^gab ^gali Mg, T , . 



' eal ' gal ' B ' 



'},(5) 



the sky, absolute B-band magnitude, and spectral type of 
the ith galaxy, respectively. 



3 METHOD 

3.1 Computing the magnification 

When computing the magnification of an individual SN la, 
we take many galaxies along the line of sight into account. 
Since the thickness of the lenses is small compared to the 
distances between them, we can use the thin screen approxi- 
mation. In this approximation the lenses are treated as two- 
dimensional surfaces perpendicular to the line joining the 
source and the observer. The surface mass density of a lens, 
E(^), is obtained by projecting the mass density, p(r), of 
the lens onto a plane (the lens plane). 



(6) 



where ^ is a vector in the lens plane and y is a coordinate 
along the line of sight. In this step of the calculations, the 
density profile of the dark matter haloes enters. 

In the weak lensing approximation the magnification 
factor, /i, is obtained by summing the convergence due to 
all lenses along the line of sight, 



^ - ^ + 2^ Kgal- 



(7) 



The convergence of a single galaxy halo, 

Kgal = S/Ec, (8) 

is the surface density scaled by the critical surface density, 
Ec — Ds/{4:TrGDiDis), where Ds is the angular diameter 



distance between the observer and the source, Di is the dis- 
tance between the observer and the ith lens, and Dis is the 
distance between the zth lens and the source. 

For most SNe la we expect the effects of gravitational 
lensing to be small. In the case of small magnifications the 
weak lensing approximation is appropriate. Since our sam- 
ple contains SNe la at very high redshift the weak lensing 
approximation might not be valid for all SNe la. We there- 
fore recursively trace a light ray through all lens planes to 
find the position of the SN la in the source plane and its 
magnification. 

The position of the SN la in the jth lens plane can be 
found from 



£2 



1=1 



-gal 



(^0 



(9) 



where dgai is the deflection angle due to the ith galaxy halo. 
Dij represents the distance between the ith and the jth lens 
plane. For spherically symmetric lenses, the deflection angle, 
which depends on the surface mass density, is given by 



(0 



4G, 



2n / m')^'d^'i, 



(10) 



where C = 1^1 ^ is a unit vector in the direction of ^. 

The magnification matrix can also be computed recur- 
sively by 



I 



E 



(11) 



where I is the unit matrix and Ai = I. In the above 
equation f3ij represents the distance combination jSij — 
DijDs/{DjDis), where Ds and Dis are the distances between 
the observer and the source and between the ith lens plane 
and the source, respectively. The matrix 



U 



Kgal + 7gal,l 
7gal,2 



7gal,2 



(12) 



depends on the convergence, Kgai, and the shear components, 
7gai,i and 7gai,2, of the lens. The magnification factor of the 
SN la is given by 

tJ-^-, — 7^ . (13) 

detAiVg,,+i ^ ' 

i.e. the inverse of the determinant of the magniflcation ma- 
trix corresponding to the source plane (j = TVgai + 1). 

In order to compute the magnification of a SN la rela- 
tive to a homogeneous universe, the magniflcation obtained 
via equation (|13|l needs to be normalised using the average 
magnification factor, {fi}. In terms of logarithmic magni- 
tudes the magniflcation of a SN la is given by. 



Amiens = -2.5 logjo + 2.5 log^o {^) . 



(14) 



To compute {fi} we use simulated lines of sight. For the sim- 
ulated lines of sight we use the weak lensing approximation, 
i.e. equation See section [5TT] for a discussion of this ap- 
proach. 

3.2 Halo model 

To model dark matter haloes we use singular isothermal 
spheres (SISs). The SIS model is characterised by a single 
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Table 1. GOODS SN la data. 



Name 


Sample" 




OSN (J2000) 


<5sN (J2000) 


AmsN (mag) 




HST04Kur 


Silver 


0.36 




32'" 


36! 03 


-27° 


51' 


17"66 


-0.40 


0.39 


HST04YOW 


Gold, Union 


0.457 


12^ 


36"^ 


34! 33 


+62° 


12' 


12^'95 


-0.01 


0.32 


HST04Haw'' 


Silver 


0.49 


12'' 


35"" 


41! 16 


+62° 


11' 


37'.' 19 


0.11 


0.24 


SN2002hr 


Silver, Union 


0.53 


03'' 


32"" 


22! 57 


-27° 


41' 


52'.' 20 


0.47 


0.27 


HST05Zwi 


Silver, Union 


0.52 


03'' 


32"" 


45! 65 


-27° 


44' 


24^' 30 


-0.54 


0.37 


SN2003be 


Gold 


0.64 


12'' 


36"" 


25! 97 


+62° 


06' 


55" 60 


-0.10 


0.25 


HSTOSDic'' 


Silver 


0.64 


12'' 


35"" 


49!61 


+62° 


10' 


11('96 


-0.23 


0.18 


SN2003bd 


Gold 


0.67 


12'' 


37m 


25! 06 


+62° 


13' 


17"50 


-0.05 


0.24 


HST04Rak 


Gold, Union 


0.739 


03'' 


32"" 


18! 15 


-27° 


44' 


10"55 


-0.12 


0.22 


SN2002kd 


Gold, Union 


0.74 


03'' 


32"" 


22! 34 


-27° 


44' 


26('90 


-0.34 


0.19 


HSTOSSpo 


Gold, Union 


0.839 


12'' 


37m 


06! 53 


+62° 


15' 


11^70 


-0.39 


0.20 


SN2003eq 


Gold, Union 


0.85 


12'' 


37m 


48! 34 


+62° 


13' 


35('30 


-0.17 


0.21 


HST04Man 


Gold, Union 


0.854 


12'' 


36"" 


34!81 


+62° 


15' 


49^06 


0.07 


0.29 


SN2003eb 


Gold 


0.90 


12'' 


37m 


15! 18 


+62° 


13' 


34" 60 


-0.39 


0.25 


SN2003es 


Gold 


0.954 


12'' 


36"" 


55! 39 


+62° 


13' 


11^'90 


0.12 


0.27 


HST04Tha 


Gold 


0.954 


12'' 


36"" 


55! 17 


+62° 


13' 


04'.'05 


-0.33 


0.27 


HST04Omb 


Gold, Union 


0.975 


03'' 


32"" 


25! 34 


-27° 


45' 


03'.'01 


-0.03 


0.26 


HST04Eag 


Gold, Union 


1.019 


12'' 


37"' 


20! 75 


+62° 


13' 


41^'50 


0.16 


0.19 


HSTOSFer 


Gold, Union 


1.020 


12'' 


36™ 


25! 10 


+62° 


15' 


23'.' 84 


-0.37 


0.27 


HSTOSStr 


Gold, Union 


1.027 


12'' 


36™ 


20!63 


+62° 


10' 


50'.' 58 


0.43 


0.19 


HSTOSGab 


Gold, Union 


1.12 


12'' 


36™ 


13!83 


+62° 


12' 


07" 56 


0.05 


0.29 


SN2002ki 


Gold, Union 


1.14 


12'' 


37m 


28! 35 


+62° 


20' 


40^' 00 


0.05 


0.29 


HST04Gre 


Gold, Union 


1.14 


Uo 


Qom 




97-0 


A f<' 
40 


00. 0\J 


-0.22 


0.31 


HSTOSRed 


Silver, Union 


1.19 


12'' 


37m 


01! 70 


+62° 


12' 


23'.' 98 


-1.14 


0.39 


HSTOSKoe'' 


Gold 


1.23 


12'' 


36™ 


22!92 


+62° 


18' 


23'.' 20 


0.30 


0.23 


HSTOSLan 


Gold, Union 


1.235 


12'' 


36™ 


56! 72 


+62° 


12' 


53'.' 33 


0.11 


0.20 


SN2003az 


Silver, Union 


1.26 


12'' 


37™ 


19!67 


+62° 


18' 


37('50 


-0.31 


0.25 


SN2002fw 


Gold, Union 


1.30 


03'' 


32™ 


37!52 


-27° 


46' 


46^60 


0.04 


0.20 


SN2002hp 


Gold, Union 


1.30 


03'' 


32™ 


24! 79 


-27° 


46' 


17"80 


-0.52 


0.30 


SN2003aj'' 


Silver 


1.31 


03'' 


32™ 


44! 33 


-27° 


55' 


06^40 


-0.05 


0.31 


SN2003dy 


Gold, Union 


1.34 


12'' 


37™ 


09! 16 


+62° 


11' 


29^' 00 


-0.17 


0.31 


HST04Mcg 


Gold, Union 


1.357 


03'' 


32™ 


10!02 


-27° 


49' 


49^' 98 


0.07 


0.25 


HST04Sas'= 


Gold, Union 


1.39 


12'' 


36™ 


54! 11 


+62° 


08' 


22^76 


-0.29 


0.19 


SN2002fx 


Silver 


1.40 


03'' 


32™ 


06! 80 


-27° 


44' 


34^40 


0.06 


0.81 


SN2003ak'' 


Silver, Union 


1.55 


03'' 


32™ 


46! 90 


-27° 


54' 


49^' 30 


-0.43 


0.32 


SN1997fr 


Gold 


1.755 


12'' 


36™ 


44! 10 


+62° 


12' 


44^69 


-0.48 


0.35 



°' The terms Gold and Silver indicate that the SN la has been classified as gold or silver by Riess et al.l 
l|2003). The term Union indicates that the SN la belongs to the compilation of lKowalski et al ] 120081 1 
^ Located too close to the boundary of the field. 
Located too close to a masked region. 



parameter: the velocity dispersion, a. Since the SIS model 
is axially symmetric, the surface density, E(^) = cr^/(2G^), 
is a function of the physical distance between the source 
and the lens in the source plane, ^. Given the distance to 
the lens, Di, and the angular separation between the source 
and the lens, 9, the physical distance in the lens plane can be 
computed by the formula ^ = Di6. For a SIS the convergence 
is given by 



Kgai(a;) — — , 

where x = ^/^o and 
2 DiDis 

4o = 47r(j 



The shear components of a SIS are 



7gal,l(x) = 



- Xi 



2x3 



(15) 



(16) 



(17) 



and 

7gal,2(x) 



X1X2 



(18) 



where xi and X2 are the two components of the vector x = 

To investigate the relationship between galaxy luminos- 
ity and velocity dispersion, we use the scaling law 



(19) 



which has two parameters: ct, and 77. L» is a characteristic 
luminosity, which we take to be 10^'^ Lq in the B-band. 
In terms of absolute B-band magnitudes, which we will work 
with, the scaling relation becomes 

-i7(A/B-i\/* )/2.5 



(20) 



where Mg = —19.52 + 51ogjQ h. Our estimates of Mb for 
the galaxies depends on h in the same way as Mg, which 
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means that h (the Hubble constant is related to this pa- 
rameter through Ho = lOO/i km s~^ Mpc~^) cancel out 
from the calculations. For GOODS-S and GOODS-N, the 
average absolute _B-band magnitude of the galaxies are 

(Mb) = -19.46 + 51ogio/i and (Mb) = -19.13 + 5 log^o ft, 
respectively. The average values of Mb are thus close to the 
value of Mg. 

Our model describes a scaling relation between lumi- 
nosity and velocity dispersion with two parameters, which 
we will refer to as 

Ahalo = W*,il}- (21) 



3.3 X -statistic 

In order to find the best-fitting halo model, with parameters 
Ahaio, we minimise the following x^-statistic: 



JVSN 

E 



Amgp^ — Amiens (-^SN, -^losj -^halo) 



(22) 



where the superscripts refer to the ith SN la. The predicted 
magnification of a SN la, Amiens(AsN, Aios; Ahaio), depends 
on the halo model Ahaio and the variables Asn describing the 
SN la position and the variables Aios describing the galaxies 
along the line of sight. 



4 RESULTS 

Figure [1] shows the result of fitting cr, and rj to the Gold 
sample. The best-fitting model, with Xmin ~ 23.71 for 22 
degrees of freedom corresponding to a reduced of 1.08 
is indicated by the circle. The solid contours show 68.3, 95, 
and 99 per cent confidence levels corresponding to x^ ~ 
Xmin + Ax^ with Ax^ = 2.30, 5.99, and 9.21, respectively. 
The best-fitting values are: cr* = 136 km s-i and r? = 0.27. 
We find CT* < 190 km s"'^ at the 95 per cent confidence 
level. Since k oc ct*, a vanishing value of cr, corresponds to 
no gravitational lensing in our model. At the 95 per cent 
confidence level, the results are therefore consistent with no 
gravitational lensing. 

If the lensing hypothesis is correct, residuals and 
magnifications should be correlated. Figure [5] shows a 
magnification-residual diagram computed for the best- 
fitting halo model for the Gold sample. We expect the points 
to scatter around the solid line, if they are correlated. The 
linear correlation coefficient is rcorr = 0. 3 6 for this sam- 
ple, which is similar to what I Jonsson et al.1 1I2OO7I) found as- 
suming a halo model base d on TuUy-Fisher iTullv fc Fisheij 
[l973) and Faber-Jackson (|Faber fc Jacksonlll976l l relations. 
For simulated data sets corresponding to no gravitational 
lensing we find r^orr ^ 0.36 for the best-fitting model for 18 
per cent of the samples. 

Our estimates of cr, and rj can be translated into a re- 
lation between luminosity and mass. The (aperture) mass 
enclosed within a radius R is for a SIS given hy M{r ^ R) = 



2cr R/G. For the scaling model described by equation (fT9 
the aperture mass can consequently be written 



d 



d 

I 




50 100 150 
[km s~^] 



200 250 



/ r \ 2)7 

M(r ^R)= Af, i^—J , 



(23) 



Figure 1. Confidence level contours (68.3, 95, and 99 per cent) 
in the (cr, , r))-plane. Circle and solid contours indicate best-fitting 
values (cr, = 136 km s~^ and r] = 0.27) and confidence level 
contours obtained for the Gold sample. Star and dotted contours 
correspond to best-fitting values (cr, = 73 km s~^ and rj = —0.18) 
and confidence level contours obtained when adding the Silver 
sample to the Gold sample. 

where M, = 2a^R/G is the meiss of a L, galaxy halo. For the 
Gold sample we find M, = 1.3 x IO^^/i'^Mq (Af, < 2.5 x 
10^^ Mq at the 95 per cent confidence level) assuming 
7? = 150/i"^ kpc. 

In order to increase the number of SNe la we could in - 
clude also SNe la classified as silver by iRiess et al.l (|2007t ) 
in our fit. The total number of silver SNe la in the GOODS 
fields is 10, but 4 of these are located to close to the bound- 
ary of the field. The residual of HST05Red is Amiens < — 1 
mag. Since there are no massive galaxies along the line of 
sight, such a large magnification is very unlikely to be caused 
by lensing. Consequently, we do not include HST05Red in 
the analysis. Our Silver sample therefore contains 5 objects. 
The best-fitting values to the joint Silver and Gold sample is 
indicated by the star in Fig.[T] The 68.3, 95, and 99 per cent 
confidence levels are represented by the dotted contours. In- 
cluding the 5 objects in the Silver sample have a very small 
effect on the confidence level contours, but leads to shifted 
best-fitting values (cr, — 73 km s~^ and r) = —0.18). 



5 SYSTEMATIC EFFECTS 

In this section we examine different systematic uncertainties 
which could potentially affect or results. Since light curve 
fitting is a particularly important source of systematic un- 
certainty, we defer the discussion of it to section [S] 



5.1 Magnification factor normalisation 

To compute the magnification normalisation, {^), we use 
simulated lines of sight. The simulated lines of sight are ran- 



6 J. Jonsson et al. 



d 



S5B 

cd 



d 

I 




-0.4 



lens 



Figure 2. Magnification-residual diagram for the best- fitting 
halo model to the Gold sample. The solid and dashed line have 
slope unity and zero, respectively. 



domly located in the image plane, but should be randomly 
located in the source plane. This leads to an over estima- 
tion of (/i) , because magnified lines of sight subtend a larger 
fraction of the image plane than de-magnified ones. To com- 
pensate for this a weighting procedure can be used. For halo 
models predicting large magnification factors, the weighting 
procedure rejects most of the lines of sight. The convergence 
of (fj.) is consequently very slow for these models. Therefore 
we use the weak lensing approximation, without any weight- 
ing, instead. We have checked that this approach is reliable. 

Figure[3]shows (fi) as a function of redshift for GOODS- 
S (solid curve) and GOODS-N (dashed curve) computed for 
the best-fitting model to the Gold sample (a* = 136 km s^^ 
and ri = 0.27). For the two fields the normalisation clearly 
differs. The normalisation depends on the magnitude limit 
of the survey and cosmic variance. For this model the weak 
lensing approximation leads to under estimated values of 
(fi) by less than 1 and 2 per cent at ^sn ~ 1 and ^sn ~ 1.5, 
respectively. The error from the approximation increases for 
halo models predicting larger magnifications. For example, 
the errors for a model with cr* = 200 km s~^ and rj = 0.27 
are 3 and 6 per cent for ^sn = 1 and zsn ~ 1-5, respectively. 

To compute {fi} for a particular model and a partic- 
ular redshift, 1000 simulated lines of sight are used. Since 
increasing this number by a factor of 2 has negligible impact 
on the results, we conclude that 1000 simulated lines of sight 
are sufficient for estimating (/i). 
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Figure 3. Average magnification factor, (^), computed for sim- 
ulated lines of sight. Solid and dashed curves correspond to 
GOODS— S and GOODS— N, respectively. The curves correspond 
to a halo model with parameters tr* = 136 km s~^ and r\ = 0.27. 



mit this approach. Instead we have to assume a cosmological 
model in our calculations. Since we restrict our investigation 
to a fiat universe with a cosmological constant, the cosmo- 
logical model can be simply characterised by Qm- Changing 
the value of Am firom 0.3, which is what we usually assume, 
to 0.25 and 0.35 results in a change in the best-fitting value 
of (T* from 136 to 152 and 129 km s~^, respectively. The cor- 
responding changes in the best-fitting values of rj are negli- 
gible. The cosmological model thus has a rather weak effect 
on the results. For the Gold sample, the effect of the cosmo- 
logical model when considering realistic deviations from the 
concordance model, is certainly smaller than the statistical 
uncertainty. 



5.3 Cut-off radius 

When deciding which galaxies to include in the analysis, we 
use an angular cut-off radius 9c. If a circle with this radius 
centred at the SN la position crosses the edge of the field or a 
masked region, the SN la is excluded from the analysis. Since 
the SIS profile is divergent, we also use a cut-off radius, ^max, 
in the lens plane. At low redshifts ^max has no effect since 
its corresponding angular scale is larger than 6c. We use 
Cmax — 150h~^ kpc. Increasing the value of ^max to 250/i~^ 
or 350/i~^ kpc has a negligible effect on the results. 



5.2 Cosmology 

Cosmological parameters enters our calculations not only 
through the angular diameter distances, but also, more im- 
portantly, through the Hubble diagram residuals. Ideally the 
residuals would be computed using the average SN la bright- 
ness, but the size of the data set we use is too small to per- 



5.4 Magnitude limit 

Our galaxy catalogues have a limiting magnitude of 25 mag 
in the i^s-band. We assume that the more luminous a galaxy 
is, the more massive and more important it is for the grav- 
itational lensing calculations. To test this assumption we 
can compare results obtained using different magnitude lim- 
its. Table [2] shows the best-fitting values of a, and rj for 
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Table 2. Effect of changing the hmiting magnitude. The table 
shows the best-fitting values of cr* and r] obtained for different 
magnitude limits in the /Cs'band. The number of early (A^carly) 
and late (Ni^^^) type galaxies in the GOODS fields brighter than 
the limit are also given in the table. 



-f^S (mag) 


cr« (km s ■'-) 


V 


N 1 

'early 




23.5 


152 


0.21 


634 


2296 


24.0 


149 


0.24 


752 


3638 


24.5 


136 


0.30 


825 


5659 


25.0 


136 


0.27 


883 


8502 



magnitude limits in the range 23.5 Ks ^ 25. The table 
also shows the number of early and late type galaxies in 
the GOODS fields brighter than the magnitude limit. The 
value of rj increases slightly when the limiting magnitude in- 
creases, while a, decreases. The gravitational lensing effect 
from galaxies too faint to be included in the calculation ap- 
pears to be compensated for by increased values of a, . Since 
the best-fitting values of the halo parameters approach the 
ones obtained for Ks = 25 mag when the limiting magni- 
tude increases and the scatter in the best-fitting values is 
much smaller than the statistical uncertainty, we conclude 
that this magnitude limit is sufficient for our study. From 
the table it is clear the number of early type galaxies is less 
affected by the magnitude limit than the number of late type 
galaxies. 



5.5 Photometric redshift uncertainties 

The redshift of the galaxies is a particularly important pa- 
rameter in the calculation of the magnification, because it 
enters the calculations not only via the angular diameter 
distances involved, but also through the estimation of the 
luminosity of the galaxies. For most of the galaxies we have 
to use photometric redshifts with considerable uncertainty. 
The error in the SN la redshift is negligible. To investigate 
the effect of the photometric redshifts on the dark matter 
constraints, we use simulated galaxy catalogues where the 
original photometric redshifts, 2gai, have been replaced by 
redshifts, z'^^i, drawn at random from the probability distri- 
bution of the photometric redshifts. A change in the redshift 
is accompanied by a correction, AMsiz'^^i), in the absolute 
magnitude. The absolute magnitude of the galaxy with the 
new redshift is hence 



= A-/B(2gal)-HAA/s(z;al) 



(24) 



Ms(4ai) 

The results of the simulations are presented in Fig. |4] Open 
circles represents the best-fitting values of cr* and rj to 100 
simulated galaxy catalogues. For the benefit of comparison, 
confidence level contours and best-fitting values (solid circle) 
obtained for the original catalogue are also shown in the 
figure. From Fig. |4]we conclude that the uncertainty in the 
results due to the photometric redshifts is small compared 
to the statistical uncertainty. 



5.6 Shifted halo and galaxy positions 

Shifts between galaxy and halo positions could potentially 
affect or results. In a study of strong gravitational lensing 
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Figure 4. Best-fitting parameters cr« and t] computed for differ- 
ent realisations of the galaxy catalogues according to the photo- 
metric redshift probability distributions of the galaxies. Note that 
the scale is different from Fig. [Tl 



systems, IYoo et all l|2006h found the centre of light to coin- 
cide with the centre of mass. Even if galaxies and haloes are 
not shifted relative to each other, measurement errors could 
lead to a shift. We have investigated the effect of shifting all 
galaxy positions by a value drawn from a Gaussian distri- 
bution with a standard deviation of 0.5". We find the effect 
of shifted galaxy positions to be negligible. 

5.7 Ellipticity 

So far we have assumed haloes to be axially symmet- 
ric. In order to test for the effect of deviations from 
axial symmetry, we make use of a generalisation of 
the SIS called the singular isot h ermal ellipsoid (SIE, 
iKormann. Schneider fc BartelmannI 1 19941 '). For the SIE 
model the convergence. 



%aU 



2x 



\J cos^ + sin^ <jf> 



(25) 



is a function not only of x — ^/^o, with given by equa- 
tion (|16|) . but also of the angle, (p, between the position of 
the ray in the lens plane and the minor axis of the lens. 
In addition to a, the SIE model depends on the minor 
to major axis ratio, /. The shear components of the SIE 
model are given by 7gai,i(a;,0) = — Kgai(a;, 0) cos(2(^) and 
7gai,2(a::,<?!') = -Kgai(a;, <^) sin(2<;/>). 

We use simulations to explore the effects of el- 
liptical haloes. For simplicity we assume all values 
of the inclination of the galaxies, ip, in the range 
< V ^ 27 r to be equally plausible. Numerical 
simulations (e.g., iGustafsson. Fairbairn fc Sommer-LarserJ 
I2OO6I) implies that dark matter haloes are triax- 



ial. 

c/a 



Gustafsson. Fairbairn fc Sommer-LarsenI (|2006l ) found 



0.73±0.11, where c and a are the smallest and largest 
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Figure 5. Best-fitting parameters cr* and r) to the SIE model 
computed for different realisations of the galaxy catalogues with 
simulated values of / and ip. Note that the scale is different from 

Fig.m 

of the principal axis of a triaxial halo. We make the conser- 
vative assumption that we are viewing all halos along a line 
perpendicular to the plane spanned by a and c, hence max- 
imising the effect of the ellipticity. In our simulations, the 
axis ratio is therefore given by a Gaussian distribution with 
mean 0.73 and standard deviation 0.11 truncated at / = 1. 
For this choice of parameters, the convergence of the SIE 
is on average 7 per cent smaller than a SIS with the same 
velocity dispersion. 

The best-fitting parameters (indicated by open circles) 
to 100 simulated data sets are shown in Fig. [S] The points 
corresponding to the simulated SIE models are scattered 
around the solid circle, which indicates the best-fitting pa- 
rameters to the SIS model. Since the scatter is small com- 
pared to the statistical uncertainty, we conclude that the 
ellipticity of the galaxies can safely be neglected. 

Although the ellipticity of the haloes is not a concern 
for the current data set, it could certainly be for larger data 
sets with smaller statistical uncertainty. In that case the 
observations of galaxy ellipticities, which are used for shear 
measurements, might be used to improve the modeling of 
the galaxies for the purpose of convergence measurements. 



the host galaxy dust corrections are based on the observed 
colours of the SNe la, they will also to some extent cor- 
rect for dust in intervening galaxies and intergalactic dust. 
However, since we do not expect this correction to be exact, 
SNe la with foreground galaxies close to the line of sight 
may suffer from systematic effects from the dust extinction 
in the intervening galaxies. This could lead to a correlation 
between the residual extinction and magnification. 

Most SNe la in the Gold sample have an impact param- 
eter of the closest galaxy along the line of sight th at is larger 
than the disk length scale (< 5^~^ kp c . iFreemanl 1 1 9 70, ) . We 
therefore expect the extinction from dust in the disks of in- 
tervening galaxies to have a small effect for most SNe in 
the current sample. For SN2003es and HSTOSLan, the clos- 
est impact parameter is 5/i~^ and 7h~^ kpc, respectively. 
Since these values are comparable to the disk length scale, 
SN2003es and HSTOSLan could be affected by disk dust ex- 
tinction. Both SNe la are ~ 0.1 mag dimmer than expected 
from the concordance model. According to the best-fitting 
halo model SN2003es and HSTOSLan should be magnified by 
—0.10 and —0.15 mag, respectively. For SN2003es the closest 
galaxy is responsible for roughly 2/3 of the magnification. 
The contribution from the closest galaxy to the magnifica- 
tion of HSTOSLan is negligible. For these two SNe la there 
is hence a large discrepancy between their measured bright- 
ness and predicted magnification. One possible explanation 
of this discrepancy is that the reddening occurs at a redshift 
different from that of the host galaxy. If these two SNe la 
are removed from the Gold sample, the best-fitting values 
of a* and r) become 154 km s~^ and —0.06, respectively. 

Since 30 < a. < 220 km s"^ at the 99 per cent 
confidence level, this reduced sample is strongly support- 
ing the gravitational lensing hypothesis. The quality of the 
fit is significantly improved when removing SN2003es and 
HSTOSLan from the Gold sample, more so than for the re- 
moval of any other two SNe in the sample. 

iMenard et al.l (|2009l) recently reported the detection of 
dust in the halos of galaxies, at distances from 20 kpc to sev- 
eral Mpc from galaxy centres, based on cross-correlating the 
brightness of high redshift quasars with the position of fore- 
groun d galaxies (see also Ostman. Goobar fc Mortsellll2006l . 
boOg ) . The dust extinction was found to have a very similar 
dependence on the impact parameter as the lensing mag- 
nification. However, since the dust extinction is an order of 
magnitude smaller than the magnification, the effect for SNe 
with large impact parameters will be negligible in this study. 



5.8 Extinction by dust 

Another potential systema tic uncertainty is dimming of 
SNe l a by dust extinction (iGoobar. Bergstrom fc Mortselll 



l2002l : iNordin. Goobar fc Jonssod 120081 ). Such extinction 
may take place in the host galaxy, in the Milky 
Way, in intervening gal a xies and in intergalactic dust 
iMortsell fc Goobad [20031 : lOstman fc Mortselll |2005| ). The 
distance moduli of the SNe la included in the analysis are 
corrected for dust extinction in the host galaxy based on the 
colours of the SN la and for dust in the Milky Way. Since 



The implications for SN la cosmology due to this de- 
tectio n was investigated in JMcnard, Kilb inger fc ScrantonI 
l|2009l ). where it was found that the effects of the intergalac- 
tic dust may not be properly accounted for by the SN la 
colour corrections, leading to biased cosmological parame- 
ters. Such a bias could be relevant for the study presented 
here, since it would affect the Hub ble diagram residuals. 
iMenard. Kilbinger fc ScrantonI l|2009t ) found the bias in the 
cosmological parameters due to dust extinction to be a few 
percent. According to section [5^ shifts in Q.m of that order 
give only small changes in the results. 
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6 LIGHT CURVE FITTING 

Light curve fitting is an important ingredient in supernova 
cosmology. Since light curve fitting affects the Hubble dia- 
gram residuals, it is also a potential systematic uncertainty 
for the analysis presented here. I n the previous ana lysis we 
have used distance moduli from iRiess et al. I (120071). which 



were obtained using the MLCS2k2 (iJha. Riess fc Kirshner 
l2007t) li ght curve fitting package. Recentlv iKowalski et al 



'2008) presented a compilation of SNe la, which includes the 
GOODS SNe la, called the Union data set. For the Union 
data set a differen t light curve fitting package, named SALT 
I Guv et al.l l2005t ) , was used. In the following we compare 
results on dark matter halo propertie s obtained using dif - 
ferent light curve fitting packages [see iKessler et all (l2009l l 
for a thorough study of the systematic differences between 
them] . 

Since iRiess et al.1 (120071 ) and lKowalski et all l|2008l) use 

different selection criteria, different GOODS SNe la are 
included in the Gold and Union samples. For example, 
SN2003b e SN2003bd, SN2003e s, and HST04Tha were re- 
jected bv iKowalski et al.l l|2008l ), because their light curves 
had i nsufficient early data, but included by I Riess et al.l 
(|2007t l in the Gold sample. A column indicating which 
GOODS SNe la belong to which sample is included In ta- 
bled IIST04Sas is located too close to a masked region and 
is for that reason excluded from the comparison below. 

Figure [6] s hows a scatter p l ot of SN la Hubble diagram 
residuals fr o m IKowalski et al.l (|2008l ) versus residuals from 
iRiess et al.] l|2007l ') for the 76 SNe la common to the Gold 
and Union data sets. We expect the points to scatter around 
the solid line (with slope unity), if there i s no systematic 
offset between the two sets of residuals. The IKowalski et al.l 
(l2008h residuals a re on average 0.013 mag dimmer than the 
iRiess et al.l l|2007l ) residuals. The scatter around the solid 
line, which is indicated by the dashed lines, is 0.22 mag. For 
the subset of SNe la located in the GOODS fields, indicated 
by filled circles , ther e appears to be a larger offset. The 
IKowalski et al.l l|2008h residuals are on average 0.055 mag 
dimmer for this subset of SNe la. Since the scatter is quite 
large, 0.20 mag, this offset might not be significant. 

Figure [7] shows a comparison betw een r esults obtained 
using residuals from lRiess et al.l ((2003) and lKowalski et al.l 
l|2008l) for the 18 GOODS SNe la in common to both data 
sets. The circle and the star corres pond to the best-fi tting 
values obtained for r esiduals from iRiess et al.l (|2007l ) and 
IKowalski et al.l (|2008| ). respectively. Solid and dotted con- 
tours show 68.3, 95, and 99 per cen t confidence level con- 
tours corresponding t o the data from lRiess et al.l (|2007t ) and 
IKowalski et al.l (|2008l ). The different residuals lead to some- 
what different results. 

Figure [S] shows a comparison between the results ob- 
tained for the Gold (circle and solid contours) and the Union 
(star and dotted contours) samples. The best-fitting value 
for the Union sample is a* = 3 km and rj = —0.64. The 
results obtained using different light curve fitting packages 
clearly disagree at the 68.3 per cent confidence level. 

A dding an offset of —0. 055 mag, as suggested by Fig. |6l 
to the IKowalski et al] (|2008l ) residuals has little effect on the 
best-fitting values of the halo parameters, but improves the 
agreement between the Gold and Union results at the 68.3 
per cent confidence level. 
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Figure 6. Sca t ter pl ot of SN la Hubble dia gram residuals from 
IKowalski et al.l ||2008|) versus residuals from iRiess et al. 
Circles represent the 76 SNe la common to the Gold and Union 
data sets. Filled circles correspond to the subset of 19 SNe la 
(including HST04Sas) located in the GOODS fields. The dashed 
lines outlines the scatter (0.22 mag) between the data sets. 
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Figure 7. Comparison of results obtained for 18 SNe la with 
different residuals. The circle and the contours show the best- 
fitting values and the 68.3, 95 , and 9 9 per cent confidence level 
for residuals from I Riess et al. I l !2007l) . The star and the dotted 
contours show the best-fitting values a nd the 68.3, 95, and 99 per 
cent confidence level for residuals from IKowalski et al 
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Figure 8. Comparison between the results obtained from the 
Gold and the Union sample. The circle and the contours show the 
best-fitting values and the 68.3, 95, and 99 per cent confidence 
level for the Gold sample. The star and the dotted contours show 
the best-fitting values and the 68.3, 95, and 99 per cent confidence 
level for the Union sample. 



The Gold and Union samples consist of 24 and 21 
SNe la, respectively. Since they have only 18 SNe la in com- 
mon, the Gold sample contain 6 objects which are not in 
the Union sample, while 3 objects are unique to the Union 
sample. Since these data sets are very small we would like to 
investigate the effects of individual SNe la in more detail. In 
order to do this we study the difference in for the best- 
fitting model to the Gold sample and the no lensing, or null, 
hypothesis. Ax' = x\Kiio) " x\Kfio)- K Ax' < the 
best-fitting model is favoured, whereas the null hypothesis 
is preferred by the data if Ax' > 0. Since Ax' is a sum, 
we can study the contribution. Ax', to it from individual 
SNe la. The diagram in Fig. |9] shows Ax' for the Union and 
Gold sample in the top and bottom panel, respectively. For 
both samples Ax' was computed using the best-fitting val- 
ues of the halo parameters to the Gold sample. The gaps in 
the diagram represents the SNe la that were not included in 
either the Union or Gold sample. This diagram reveals some 
important differences between the two data sets. For the 
Gold sample Ax' = -2.81 with HST05Str, HST05Fer, and 
SN1997ff giving large negative contributions. HSTOSLan and 
HST04Eag, on the other hand, give large positive contribu- 
tions. For the Union sample only HSTOSStr and HST05Fer 
give large negative contributions, because SN1997ff is not 
included in this sample. HST05Lan gives a large positive 
contribution to Ax' also for the Union sample. It is clear 
from the diagram that SN1997ff is responsible for much of 
the difference in Ax' between the two samples. The results 
thus depends on the contributions from individual SNe la, 
to a relatively large extent. 

Figure [To] shows a scatter plot of Ax' versus Amiens for 
the Gold sample. Since there is no sign of any correlation 
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Figure 9. Contributions from individual SNe la to Ax^ = 
x'(-^hllo) - x'C-^Sai" ) tlif^ Union (top panel) and Gold (bot- 
tom panel) samples. The best-fitting model, ^JJ'jq, to the Gold 
sample was used to compute the differences, Ax^. 



between Ax' and Amiens, SNe la contribute both positively 
and negatively to Ax' regardless of whether they are mag- 
nified or de-magnified. 

In section [FiSl we discussed the possibility of HSTOSLan 
being affected by dust extinction. Extinction by dust could 
be the reason why HSTOSLan opposes the gravitational lens- 
ing hypothesis. Excluding this SN la from the Union data set 
leads to best-fitting values cr, = 84 km s~^ and r) — —0.15, 
which are in slightly better agreement with the gravitational 
lensing hypothesis. 



7 SUMMARY AND DISCUSSION 

Distant SNe la at high redshift might be affected by gravita- 
tional lensing due to matter in the foreground. Gravitational 
lensing magnification or de-magnification should correlate 
with Hubble diagram res iduals. By exploiting this correla- 
tion (| Jonsson et al.|[2007l ). we measure the relation between 
galaxy luminosity and the velocity dispersion of dark matter 
haloes using data from the GOODS. 

For a scaling law, a = a,{L/Lf)^ , describing the re- 
lationship between galaxy absolute B-band luminosity, L, 
and dark matter velocity dispersion, we find the best- 
fitting parameters to be cr, — 136 km s^^ and rj = 0.27 
for L, = lO^°/i"'L0. These results agree well with the 



results from galaxy -galaxy lensing studies. I n a galaxy- 
galaxy lensing study iKleinheinrich et al.l l|2006h found ct, = 
I32I24 km s~^ and rj = 0.37 ± 0.15 for the same scal- 
ing la w and halo model as we use here. iHoekstra et al.l 
(j2004h also used galaxy-galaxy lensing to investigate dark 
matter, but used the truncated iso thermal sphere model 
jBrainerd. Blandford fc Smaill Il99d ). which has an addi- 
tional parameter describing the truncation radius, s, of the 
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Figure 10. Scatter plot of contributions from individual SNe la 
to Ax^ = X^(-^haio) ~ -^^(■^hal'o) ^f™^ magnification computed 
for the best-fitting model for the Gold sample. 



halo. For a value of 77 fixed to 0.3 they found a* = 136±5±3 
km (statistical and systematic uncertainties), which 
agrees with our result, and s = 185^28 '^^^ '^P'-- 

The constraints derived from galaxy-galaxy lensing are 
stronger than the ones derived here. Considering how small 
our sample of lensed SNe la is, it is encouraging that we can 
derive any constraints at all. Our results rely upon a dif- 
ferent aspect of gravitational lensing and therefore provide 
an independent cross-check of the galaxy-galaxy lensing re- 
sults. 

We have also investigated different systematic uncer- 
tainties. The effects of the uncertainties in the cosmological 
model, photometric redshifts, and magnitude limit appears 
to be small. Light curve fitting is another systematic effect, 
which unlike the other ones cannot be ignored for the data 
set used here. Different light curve fitting packages yield dif- 
ferent distances for the same light curves. The Gold and 
Union samples, which are both small and include different 
GOODS SNe la, give rather different results. The Union 
sample is consistent with the no lensing hypothesis (cr, — 0) 
at the 68.3 per cent confidence level. The Gold sample is 
consistent with the no lensing hypothesis at the 95, but not 
at the 68.3 per cent confidence level. From Fig. [^it is clear 
that much of the differences can be attributed to individual 
SNe la. Clearly, light curve fitting is an issue which has to 
be resolved before SNe la can be used to their full potential 
as probes of the dark side of the Universe. 

The effect of removing two SNe la from the Gold sample 
was shown to be quite large in section [5^ For this reason we 
have performed a bootstrapping exercise to investigate the 
robustness of our results. The open circles in Fig. 1111 show 
the best-fitting values for 1000 simulated data sets based 
on resampling of the Gold sample. To facilitate comparison, 
the best-fitting values (indicated by the solid circle) and 
confidence level contours for the original sample are also 
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Figure 11. Result of bootstrapping exercise. Open circles cor- 
respond to the best-fitting values of cr* and rj obtained for 1000 
simulated data sets obtained by resampling of the Gold sample. 
For comparison the best-fitting values (indicated by the filled cir- 
cle) and confidence level contours for the original Gold sample 
are also shown. 

plotted in the figure. The highest density of points is not 
found near the best-fitting values of the original sample, 
but closer to ~ -0.2 and 30 < cr. < 120 km s"^ This 
is another indication that our results depends on individual 
SNe la. 

From our study it is clear that more high redshift SNe la 
are needed to obtain meaningful results. We anticipate that 
the final data from the Supernova Legacy Survey (SNLS, 
lAstier et a l. 2006) will be ideally suited for an investiga- 
tion o f the type presented here. According to lJonsson et al.l 
(|2008t ) a firm detection of the correlation between gravita- 
tional magnification and SN la brightness can be expected 
for the final SNLS data set. The SNLS produces not only 
a homogeneous set of well calibrated SNe la, but also the 
necessary deep multi colour observations of the foreground 
galaxies. 
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